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Abstract 

A method to- estimate the_frequency. of- transmi-s.- 
s.lon. overhaul s i t presented. Thi s. method, fs based 
on the two-parameter WeibuJT statistical dlstrib-- 
tton for component Fife.. A second, method Is 
presented, to. estimate. the.tiumber of reo-lacement com- 
poneats needed, to. support the transmission overnauL 
pattern— The second-method Is based oh renewal 
theory. Confidence statistics are applied with 
both methods to Improve the statistical estimate if 
sample behavior. A transmission example is also 
presented to illustrate the use of the methods. 
Transmission overhaul frequency and component 
replacement calculations are-included in the 
example. 

Nomenclature 

b Weibull slope 

e base of the natural logarithm 

F- probability distribution function, probability 
of failure 

F k probability of at least k failures 

f probability density function 

Ln natural logarithm 

ft life, hr 
M renewal function 

M e approximate renewal function 

N r number of replacements 

0 sample size 

R probability of survival, (1 - F> 

x Integration time- variable, hr 

2)0 number of standard deviations from the mean 

which cuts off a 10-percent population tail 

T the gamma function 

6 characteristic life, hr 

P 3 third moment of a probability density function 

o standard deviation 


Subscript su_ 

a v avecage_or_ mean 

f Weibui-1 function 

t inoex 

k index 

me approximate renewaJ-function 
n number of" components In system 

r replacement function 

s system 

10 90-percent reliability 

90 9&-percent confidence 

Introduction 

The in-flight service reliability of aircraft 
transmissions is much greater than the design relia- 
bility of their components. Transmission overhauls 
provide the difference. By monitofing the onset of 
potential transmission fatigue failures, just-in- 
time overhauls maintain the transmission economi- 
cally. • One cause for propulsion system overhauls 
Is ttie finite fatigue life- of drive system conipo- 
nents. The-two-oarameter Weibull distribution 
^crlbesjhrstattstlcs-of-drive system bearing and 

Component reliabilities and lives affect trans- 
mission maintenance costs which are significant. 
Estimates of these costs are Important in the design 
stage of a transmission. 6 The two-parameter Weibull 
distribution provides information on component reli- 
ability and life. It does not predict overhaul fre- 
quency directly. 

Two steps are required to convert component 
life statistics into overhaul frequency values. 

The first Is the transmission system life model. 

This model Is a two-parameter Weibull distribution 
for the transmission system life, '.8 The second 
step is renewal theory. 

Renewal theory Is a statistical model which 
describes the maintenance cycle. The theory consid- 
ers the ongoing sequence of: use, failure onset 
repair, and return to use. For this sequence, 
renewal theory predicts the frequency of component 
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replacement ana the number of replacements needed 
to support the service maintenance schedule. *-'2 

ma< . h J? nf ! a ! nce J heor/ - sm 0 > ements these statistical 
methods wUft- estimates of the likelihood of the pre- 
dictions. Higner confidence levels require more 

s S ?tIItiX^?'f 0Ver 4 9reW6r r4n9e °^ P0SSible 

The purpose of the research presented, is to 
n JJ ethod ? 1 ° 9y / or " « 1 cu 1 atl ng transmission 
ire ana the number of comoonent replacements. The 
paper presents the theories and applies them to a 
simple transmission to illustrate their use- Esti- 
mates.^ drive system and component lives and 
replacement needs are essential in design. These 

P rov ^ a comparison of the re 1 at 1 ve-worth 
of different designs from a-safety and maintenance 
cost perspective. They also help-assess. the cost of 
operating a proposed drive system. 

Component Life and- Reliability 

The two-parameter Weibuil dt stributioa. i s com- 
•non ly used to descriUe-fatlgue life data, rt can 

J wtde variet ^ of ,ife Patterns. The reli- 
aoility of a component is the comolement of Its 
pro&abi 1 i ty of fai lure . 

Jn statistics, reliability is a double nega- 
tive. Reliability or the act of surviving is the 
state of not having failed. Statistics count 
direct events such as the act of failing. A oart 
can fail only once. It survives for its entire 

lilt; J nus prooao * ’ * ty of failure is a direct 
statistic. The probability of failure for the two- 
parameter Weibuil distribution is: 


figure I Is a plot of the ratio of the •narar- 
teristic life to the-90-oercent <-ellabil'ty i<fe i3 

a f “" ct i 1 ibu " s!coe - Substituting V 

of Eq. (4) in Eq. (3) gives: ’ 


W (J) ■ La (oXii^' 


I 


( 5) 


• 0 




(i) 


to are^nf ?K a . < ! > 5 the fo used - W manufacturers 
to present the two-parameter Wei bull ois-tti button 
characteristics of bearings..’ 3 

rsi iahin+ t ^ Eqs - the logarithm of the 

ra m ]s- proportional to the life 

of either 1 ope - TaMn 9 the logarithm 

°l father equation generates a- straight line slot 

" Fl9 ” 2 ’ The P ,ori S-s probabiiitv grapn 
for the two-parameter WeUjuU distribution. ‘ 3 P 

This graph-aids in determining the disrribu- 
l °? f P H r ! meter va,ues - for fatigue test data. 14 The 
are the resu,t ^ of a- series of- 
dentical life tests for a sample- set of hdenti^ar 
components. The f+rst failure determines the h^gh- 

t?n-« Jhi 4 ** M ! y , data point> The ne ^ failure de?er- 
m^nes the next lowest reliability data point, and 

(MTT F I he ?r e : 494 ,ife ,s the Mean T,me to Failure 
„ rt is th * sum of 4,1 times to failure 
divided. by the total number of the failures. The 

failures for a continuous probabil- 
ity d stribution is unity by definition. The <jm 

nrJ! V'!? e r t0 fai1ure ls the integral of the 

function° f Tha e i? r il tfe 4nd the P r °Pabllity density 
" ?• The T ' imits on the Integral are from ie-o 

to infinity. The average life is: 


T"*™. F 1 ,S P r °Pabt 1 1 ty of failure expressed as 
a decimal, e is the base of the natural logarithm, 
ft is the component life In million load cycles or 
hours e- is the characteristic life in million load 
p y pl ®* ® r ,t |0urs ' and P is the Weibuil slope. The 
two Weibuil parameters are 9 and b. 

The derivative of Eq. <1) with resoect to life 
is the probability density function, f: 

t 9 [ej e <2> 

If e 1 ?S ob ! b !! ,t X dens 1 ty function is a histo- 
g am of life failures for a unit population It" 
presents the scatter in the component ilves. 

The Weibuil reliability function Is often 
expressed as: 


av 


MTTF . F 4f(si)(ja. 
J 0 


( 6 ) 


tlon to rwi* b f ure distribution, the solu- 

funrtfon^r IS h??'' 44 ' nv0,ves the weil known gamma 
lu’Ji 5 t The so,utlon Is the gamma function 
multiplied by the char-acteri stic life, 9 . 


av 


MTTF 


9r (' ♦ e) 


(7) 


nfo fl 9 ^ 8 * ,s 4 P ,ot of the ratio of the average 
life to the 90-percent reliability life as a func- 
t on of the weibuil slope. The average life equals 
the characteristic life when b . 1 . 

tlon Iv standard deviation of a failure distrlbu- 


“® • ifi 


s: 


<ft - ft av ) , ’f( , i)d2 


1/2 


f 8) 
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where R is the probability of survival <1 - F) 
For a 90-percent probability of survival. R - 0.9 
g"ves^ " 4,0 So,vin 9 tor the characteristic life 


In terms of the characteristic life, the Weibuil 
tbe 9 amma function, the standard devia- 
tion of the two-parameter Weibuil distribution is- 


0 ■ Ln 


/j\- ,/b 
'.0.9/ 
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(4) 


The standard deviation of a distribution is a 
measure of the scatter of the distribution! It it 
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average life. Figure 4 is a plot of the ratio of 
the standard' deviation to tne 90-percent reliabil- 
ity life as a function of the Weibull slope. At a 
slope of b ■ 1, the Weibull distribution is the 
exponential distribution and has a large scatter. 
As the slope Increases, the scatter decrease-s- 
rapldiy. 

System Life and Rellabll .v 


Renewal Theory 


In addition to scheduling maintenance perisos, 
a service procedure musvalso est’mate tne numper 
of replacement components required. Renewal theory 
adds the renewaJ- function to the stat'st'cal tools 
for estimating repair. Ifestimwes the numoer sf 
replacements as a function of a component failure 
distribution and Its life.® - " 


One model for the life of a drive system is the 
strict series probability model. 8 This model com- 
pares a system of load, carrying gears and bearings 
to a chain of links. A chain.falls.wheh any single 
rink fa-1 Is. So too, a drive system requires repair 
when any component requires replacement or repair. 

In the strict series probability model, the relia- 
bility of a system, R s , Is the-product of the relia- 
bilities of all the components. 


Renewal theory assumes the replacement of - 
failed components-when they fail. This models an 
unending sequence of use and. repair. Aircraft drive 
system maintenance follows this pattern closely. 

The renewal function results from a sequemce-of-sta- 
tistlcally predicted, fai lures. 

Consider the maintenance sequence. In a given- 
life periods any number of failures may occuc- The 
probability- of at least one failure within a gi /en 
life from the start-ot-cperatlon is: 


n 


R, - E. R 

* 1.1 


< 10 ) 


The high speed of drive system components and- 
the scattering of loose debris warrant the strict 
series probability model. If any component fails, 
debris may be present which could damage other com- 
ponents. Therefore, a drive system requires an 
overhaul to return it to a high state of reliability 
once any element fails. 


Fj-U) . ?U> 


f f(x)dx 
J 0 


< 13) 


The probability of at least two sequential 
failures In the period is the probaoliit./ of two 
independent events. The first component must fall. 
Then a second component must begin its service life 
at this failure life and also fall. The prooabiUty 
of naving at least two failures in this period is: 


Taking the logar 1 thin of the reciprocal of i 

Eq. < 1 0 ) and using Eq. (5) for each component F,U) . f F,(a - x)f(x)dx 

yields: * j. l 



( 11 ) 


In Eq. (11), D s is the life of the entire drive 
system for the system reliability. R s . It is also 
the life of each component at the same drive system 
reliability, R«. For consistency In Eq. (11), all 
the component lives must be defined In the same 
units. The unit chosen i; hours. 


Equation (11) Is not a simple two-parameter 
Heibull relationship between system life and system 
reliability. The equation is a true two-parameter 
Welbuil distribution only when all the Weibull expo- 
nents, Dj , are equal. In general, this is not the 
case. Thus, R 5 as a function of l s when plotted 
such as in Fig. 2 may produce a curve rather than a 
straight lire. A true two-parameter Weibull distri- 
bution can be approximated quite well, however, by 
fitting the curve using a least squares method. 8 
The slope of the fitted straight line Is the drive 
system Weibull slope, b s . The life at which the 
drive system reliability equals 90 percent on the 
straight line is s. 5 iq- The Drive sy.tem two- 
parameter Weibull relationship is then: 



In Eq. (14), x is the time at wnich the first 
failure occurs. This can happen any time between 
zero and 1. At x . ), the entire life *s avail- 
able for the first failure probability. The proba- 
bilities of the second failure and the comoi nation 
event are both zero. As x Increases from zero to 
1, the probability of the first failure happening at 
time x decreases. The probability of the second 
failure Increases. At x-. a, the entire life is 
available for the second failure probability. The 
probability of the first failure is zero. The prob- 
ability of the combined event Is thus zero as well 
at x - a. The Integral defines a function for the 
probability of at least two failures in the life 
period from zero to a. 

Equation (14) repeats Indefinitely with 
increasing subscripts. The probability of having at 
least k failures In tne period from zero to a 
is: 




P k-1 ( 4 “ x x 


M5> 


In Eq. (15), Ffc-jU - <) is the probability of hav- 
ing at least k - 1 failures in the period from zero 
to % * x. The probability of at least k failures 
increases as the number of failures decrease. Also, 
the more life available for a failure, the greater 
is the chance that it .-/ill occur. 


% 


3 


The mean number of failures Is the Infinite sum 
of the probabilities of at least k failures in the 
life period, 'l. This function, MU), Is the renewal 
function. It is expresses as: 

• 

MU) - £ F k u> (16) 

k«1 

Equations (IS) an<M16> yield the number of* 
replacements, needed, to support a-maintenance sched- — 
ule. The solution Involves a- large series of 
convolution integrals-. The equations apply to any 
fa-iluc^ d-i*stributiorv. However, the solution Is. not 
easy to obtain. Figure S snows the renewa-1 function 
for a component with a two-parameter Wei bull relia- 
bility. The component life has 0 « 5000 hr and- 
b « 1.5. Tabulated solutions to the renewal, func- 
tion- for the twcr-parameter We4bu-U- disXciiwtl on-are 
available J ' *’ 4 

An aoproximation for-the renewal function from 
Ref. 10 is: 


M (U-* 
6 


a 





(17) 


The approximation accuracy increases as fc 
increases. Equation (17) is an asymptote to ttie 
exact renewal function for low-scatter distribu- 
tions. For high-scatter distributions It approxi- 
mates the renewal function closely. 


The standard deviation of the renewal function 
gives a measure of the scatter in replacement needs 
from one sample to tne next. Figure 6 is a plot of 
the renewal function standard deviation versus life 
for a component which has a two-parameter Weibull 
reliability distribution. It has a characteristic 
life of 5000 hr and a Weibull slope of 1.5. 


The approximation for the standard deviation 
of the renewal function 1 s :‘0 



number of drive systems under serv'ce J : a Hmtteo 
sample. Confidence statistics estimate hew 2 ? **er- 
ently \ small sample may benave :: -nive^sal 

population. It uses trre-snndara :ev*arcn of “he 
universal failure distribution ano pe :amcie :* 2 e 
to estimate the mean of the sample .'* 1 

For many samples of the same size, :ne mean of 
the samples has a normal dl s-tribution aocut the 
overall mean. The s tandard-de^ia tlon-of the means 
rs: 


** — $ 

where Q. is the size of the sample. 

The standard deviation of me numcer of 
replacements i$~: 


o r * lS_7me ^ 1 ) 

Arl so the total number of reel aceme r ts f or tne 
sample is: 

fir . QMU) ( 22) 

In reliability predictions, the ’ower confi- 
dence bound has much significance in aircraft 
applications. Systems that are less reliable man 
the average are important to identify for safety 
and economical reasons. The confidence distribu- 
tion estimates the mean life which will be lower 
than the mean life of a chosen percentage of all 
samples of a given sl 2 e. This life is ‘ess than 
the mean life for the entire population. For a 
90-percent confidence: 

*av,90 * ‘*av “ 2 lQ c av 

where z\q is the number of standard deviations 
below the mean which cuts off 10 percent of the 
population. For a normal distribution, 

2^0 ■ 1-282. 4 ** 4 90 percent of the normal distribu- 
tion lie above &av,90 and ^ percent Me be-low 
*av,90* 

With a 90-percent confidence that the replace- 
ments will be less, the replacement estimate for a 
component from 2 ero to life l is: 


(18) 


Nr ,90 ■ N r ♦ 2 1 Q a r 


C4>- 


where 1*3 is the third moment of the life distribu- 
tion. For the two-parameter Weibull distribution, 
the third moment is: 


a 3 m)da * ©■ 


' r (' * D 


(19) 


Since the behavior of samples differs from the 
behavior of the "ideal" distribution, conf’dence 
estimates are helpful. With the confidence esti- 
mates. one can see during the design anase the 
effects of sample s 1 ze on the life and "eclacement 
estimates. 


Figure 7 is a plot of the ratio of the third moment 
to the 90-percent reliability life as a function of 
the Weibull slope. 

Confidence Statistics 

In predicting replacement rates and maintenance 
inventories, direct theory provides mean or average 
estimates. These estimates come from the statistics 
of a universal population (that is, an infinite 
number of samples). In any real situation, the 


Example 

Mean Life 

Consider the single mesh transmission shown *n 
Fig. 8 . Assume the 90-percent reliability lives for 
the bearings and gears are given as those snown in 
Table 1. Also assume a Weibull slope of 1.2 for the 
bearings and 2.5 for the gears. It is desired to 
determine the transmission mean life with 90-percent 
confidence for a fleet of 0 * 50 aircraft. 
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From Eq. (7) Dr Fig. 3, the average Mves *ere 
determined for eacn component. C r*cm £q. > 9 ) 

Fig. 4- the standard deviations for eacn component 
*ere determined, "he results are shown in “aoie i. 


me 


M 3 560)" , A 

■ !6 :co > 3 


j Ul :oor • 1 3 S 60 > : i 
L 4d6 j 


Hie transmission system 90-percent reiiaoility 
iifs was determined based on the component lives 
using Eqs . (li) and (12). Fitting Eq. (12) to 
Eq. (11) for this da-tauyields a two-parameter 
Welbull slope of b$ * 1.57 for the transmission. 

The system 90-percent reliability life 4s 
10 * 1060 nr. From Eq. i4) or Fig. | t the trans- 
mission characteristic life is 6 « 4440 nr. From 
Eq. (7) or Fig. 3, the transmission average life 
is 2a v • 3990 hr. From Eq- (9) or Fig. 4, the 
standard deviation of-the transmission life is 
of • 2600 hr. Table 1 includes these results. 

With this data, one can estimate the overhaul 
frequency for a fleet* of similar aircraft. The mean 
life of a small sample of aircraft could be lower 
than that" for an infinite population of aircraft. 
Confidence statistics estimate the sample properties 
from the universal properties. 


! 3c 16 :90r . 5(13 560)" : 


2 ' * I 

2< ‘ 5 2CO) 3 ( 



, 1/2 

0.1038) 


( 29 > 


From ta. v 2 1 > . the standard deviation of me numce r 
of ^eolacemorits is: 



f i 
Ml 



' 30) 

with a 90-oercent confidence chat -tne replacements 
wi 1 i oe ess, the replacement estimate *or cear ; nq 
one from oe:*o to life i using Eq. (24>-4-s; 


Let us estimate the overhaul frequency for 
the fleet which use these transmissions with a 
90-percent confidence that the frequency is lower. 
From Eq . (20), the standard deviation of the mean 
life's distribution is: 


2600 
av * yso 


363 hr 


f 25) 


Using Eq. <23). the estimate of the transmission's 
mean life witn a 90-percent confidence is: 


' A av , 90 s 39 90 - 1.282(363) - 3520 hr (26) 
Number of Replacements 

The renewal function serves to estimate the 
number of components one needs to support the main- 
tenance pattern. Consider bearing number one in the 
single mesh transmission of Fig. 3 and Table 1. It 
is desired to determine the number of replacements 
required at any given time with 90-percent confi- 
dence for a sample of 50 aircraft. 


From Table 1, the average life is 
^av 2 lb 2C0 hr and the standard deviation is 
°f‘ * 13 560 hr. The renewal function from Eq. (17) 


■ ; 90 


354 * 7 ‘ 5 ) * -2821 4 -5.19j 


. 1/2 
} 


(31) 


As an e<ample, -et us use this information to 
determine the numoer of replacements '-squired for 
bearing cne if the 50 aircraft are to operate up to 
l = : 0 XO hr. Note that if **enewai tnecr/ is nor 
used, tne orooaoility of ^ai’ure from Eq. : 1 ) is 
F * 0.41 -or bearing one at i « 10 000 hr. This 
wcul. 'eao to 50(0.41) « 21 bearing required, "his 
method, holies that bearings *or 21 aircraft have 
failed ano these aircraft are no longer in service. 
The renewal theory method, however, imoMes that the 
failed Dealings are replaced and tnese' aircraft are 
Put bacK in service. For i . 10 000 hr. Eg. < 3 7 > 
estimates ttie need for 30 bearings to sueport the 
overhaul needs of the 50 aircraft, compared to 21 
using only tne Weibull failure di striouticn. 


Summary of Res u 1 ts 

A metnod to estimate the frequency of transmis- 
sion overhauls was presented. A second method was 
presented to estimate the number of replacement com- 
ponents needed to support the transmission overhaul 
pattern. Confidence statistics were applied with 
both methods to improve the statistical estimate of 
sample behavior. 


* rrroo 


(1 5 :0Q>* - 03 550)- 
2(1(3 200) 2 


V, 

16 200 * 


0.150 


(27) 


From Eq. (22), the total number of replacements from 
0 hr to a 1 i f e of % is; 


f V - 50 (rrloo - °- ,so ) ■ m - 7 - 5 

From Eq. (19) or Fig. 7, the third moment of the 
bearing life distribution is ^3 * 1 . 697 x 1 0 1 3 hr 3 . 
From Eq. (13), the standard deviation of the renewal 
function for bearing one is: 


T he metned to predict overhaul frequency is 
based on a two-parameter Weibull system life 
model, "he relationship between the system life 
model and the component life models was presented. 

In addition, formulas *or the mean and standard 
deviation of the two-parameter Weibull distribution 
were given. 

Renewal theory was presented as a ted to esti- 
mate the number of component replacements in a 
transmission. Approximation formulas were given 
for the mean ana standard deviations qf the renewal 
function. These approximations are valid for tne 
cwo-parameter Weibull distribution. Formulas for 
sample replacement rates were given in terms of the 
renewal function. 


5 


Single sided confidence theory m as presented 
for the overhaul frequency and tne :omponent 
reolacement rate esthetes. A transmission example 
was aUo 'presented to illustrate the use if the 
methods. Transmission overnaul frequency and compo- 
nent replacement calculations were included in the 
example. 
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TABLE 1. - SINGLE MESH TRANSMISSION 
PROPERTIES 


Bearing 1 


Bearing 2 
Pinion 


Bearing 3 
Bearing 4 


Transmission 1060 


i 12 




WE! BULL SLOPE, b 

FIGURE 1. - CHARACTERISTIC LIFE TO 90-PERCENT RELIA 
BILITY LIFE RATIO FOR A WEIBUU DISTRIBUTION AS A 
FUNCTION OF THE WE I BULL SLOPE. 

















FIGURE - STANDARD DEVIATION TO 90-PERCENT RELIA- 
BILITY LIFE RATIO FOR A WE I BULL DISTRIBUTION AS A 
FUNCTION OF THE WE I BULL SLOPE. 



LIFE. 1 . HR 


FIGURE 5. - RENEWAL FUNCTION FOR A TWO-PAR AHETER 
WEIBULL DISTRIBUTION WITH 0 - 5000 HOURS AND 
b = 1.5. 
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THIRD MOMENT. / (90-PERCEHT RELIABILITY LIFE) 3 , 



LIFE, t. hr 


FIGURE 6. - RENEWAL FUNCTION STANDARD DEVIATION FOR 
A TWO-PARAMETER WEIBULL DISTRIBUTION WITH 0 = 500 
HOURS AND b = 1.5. 
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WEIBULL SLOPE', b 

FIGURE 7. - THIRD MOMENT TO 90-PERCENT RELIABILITY 
LIFE RATIO FOR A WEIBULL DISTRIBUTION AS A FUNC- 
TION OF fHE' WEIBULL SLOPE. 
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